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O 



We report a recent calculation of the heavy quarkonium energy levels and decay widths in a 
quark- gluon plasma whose temperature is much smaller than the inverse radius of the bound 
state, based on a Non-Relativistic Effective Field theory framework for heavy quarkonium at 



finite temperature. Relevance for the phenomenology of the Y(1S) in heavy ion collisions is 
also discussed. 

o 

Oh 1 Introduction 

The suppression of quarkonium has been hypothesized 25 years ago [1] to represent a 
signature of the formation of a deconfined medium and has been ever since intensely 
^ investigated, both theoretically and experimentally. Here we address the problem, central 

to these studies, of the behaviour of a quarkonium bound state in a deconfined thermal 
medium. To this end we shall illustrate the Effective Field Theory (EFT) framework that 
has been recently constructed in [2-4] (see also [5-7] for an analogous EFT of QED) by 
generalizing the successful zero-temperature framework of Non-Relativistic (NR) EFTs 
for heavy quarkonia to finite temperatures. These NR EFTs exploit the hierarchy m 3> 
mv » mv 2 that characterizes any NR binary bound state, m being in this case the heavy 
H quark mass and v the relative velocity, mv is then the typical momentum transfer or inverse 

j> radius and mv 2 ~ E the binding energy. The low-lying quarkonium states, especially the 

bottomonium ground states Y(1S) and rjb, are believed to be approximately Coulombic. 
That corresponds to having mv ~ moc s 3> Aqcd arid mv 2 ~ mot 2 > Aqcd- 
In a weakly-coupled plasma, which we consider in our study, the temperature T and the 
chromoelectric screening mass jhd are larger than Aqcd and, since m 2 D ~ g 2 T 2 , T 3> mrj. 
Under these conditions one can then calculate observables relevant for the phenomenology 
of low-lying states to a large extent analytically in perturbation theory, which makes them 
extremely interesting, also in the light of the recent CMS measurements of the suppression 
of the Y family [8]. 
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In [9] the perturbative QQ static potential was computed in QCD for distances r such that 
T 3> 1/r > mi). The resulting potential surprisingly shows an imaginary part which is 
larger than the screened real part for 1/r ~ mo. This imaginary part can be traced back to 
the imaginary part of the gluon self -energy and is due to the Landau-damping phenomenon; 
it eventually leads to a thermal width for the bound state, which is in turn responsible 
for its dissociation, representing a change from the previous colour-screening paradigm. 
This change was further reinforced by the introduction of a dissociation temperature in [5, 10], 
defined as the temperature for which the imaginary part of the potentials becomes of 
the same size of its real part; parametrically it is of order ffl«P and hence smaller than 
the temperature at which screening sets in. A quantitative calculation of the dissociation 
temperature for the Y(1S) can be found in [6] and a phenomenological analysis, based on 
these imaginary parts, of bb bound states at LHC energies can be found in [11]. 
In [2] the static QQ was first studied in an EFT framework, systematically exploring the 
hierarchy of different energy scales in the problem. Many possibilities were considered, from 
temperature smaller than E to temperatures much larger than 1/r, where the results of [9] 
were recovered in a rigorous EFT derivation. Furthermore a new dissociation mechanism, 
the colour-singlet to octet decay, was identified; it is the leading one when E 3> mo- In [12] 
the relation between the proper real-time quarkonium potential and the correlator of two 
Polyakov loops, a quantity often measured on the lattice and used as input for potential 
models, was investigated. The breaking of Lorentz invariance induced by the preferred 
reference frame introduced by the medium was instead analyzed in [4] in the spin-orbit 
sector of the EFT. In the following we will report about the findings of [3], where in a specific 
range of temperatures the spectrum and width of quarkonia have been computed up to 
order ma 5 s . To this end, the specific global hierarchy that the NR and thermodynamical 
scales fulfill in the assumed range of parameters has been exploited by constructing a 
corresponding tower of EFTs. 

2 Energy scales, the EFT formalism and the results 

The aforementioned global hierarchy we assume ism > moc s 3> tzT 3> mal 3> mo, which 
implies a temperature below the dissociation temperature. We also remark that for the 
Y(1S) at the LHC it may hold that m b » 5 GeV > m b a s » 1.5 GeV > nT « 1 GeV > 
mocl 0.5 GeV > m D . 

Given this hierarchy, we now proceed to integrate out each scale in sequence. The integra- 
tion of the mass scale m yields non-relativistic QCD (NRQCD) [13] and the further integra- 
tion of the scale ma s from NRQCD gives potential non-relativistic QCD (pNRQCD) [14]. 
Since the temperature is much smaller than these two scales, it may be set to zero in the 
matching and the Lagrangians of NRQCD and pNRQCD are the same as at zero tempera- 
ture. 

Integrating out T from pNRQCD modifies it into its Hard Thermal Loop (HTL) version, 
PNRQCDhtl [2, 15], where the light degrees of freedom are described by the HTL effective 
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Lagrangian [16] and the pNRQCD potentials receive a thermal part. Finally, within this 
EFT we can compute contribution to the spectrum and width from the scales E and m D . 
Diagrams contributing to the calculation are shown in Fig. 1. 




Figure 1 : The diagrams contributing to our calculation. Single lines are colour-singlet QQ 
states, double lines are colour octets, curly lines are gluons, vertices are chromoelectric 
dipoles and the blob is the gluon self-energy. The imaginary part of the first diagram 
yields the singlet-to-octet decay mechanism, whereas the second one gives the Landau 
damping contribution to the width. 



Let us now show the final results for the thermal contribution to the spectrum and to the 
width up to order trial. We recall that for a Coulombic bound state the spectrum is at LO 
given by the Bohr levels E n = —mC 2 0L 2 / (4n 2 ) and the Bohr radius is flo = 2/ (mC F a s ). The 
vacuum contribution to the spectrum up to order m&\ can be read from [17]. 
The thermal contribution to the spectrum then reads 
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where L„ / is the QCD Bethe log [17]. The terms on the first line are the leading ones the 
power counting of the EFT and, being positive, lead to an increase in the mass of the bound 
state quadratic with the temperature. 
For what concerns the thermal width, we have 
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where I n j is a new Bethe logarithm [3]. The terms in the first two lines are the leading ones 
and are caused by singlet-to-octet decay, whereas those on the last two lines are due to 
Landau damping. The width is at leading order linear in the temperature and much smaller 
than the binding energy. This small width is certainly not in contradiction with the recent 
CMS results [8] that point to a substantial survival of the Y(1S) at the LHC. 
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